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SUMMARY 


An exploratory analysis has been developed for the case of distrib- 
uted injection of a foreign gas into a turbulent boundary layer in air 
on a flat plate. The work is divided into three parts: a derivation of 

the basic turbulent boundary-layer equations for a binary gas system; a 
derivation of modified Reynolds analogies between momentum, mass, and 
heat transfer for a binary gas system; and an evaluation of the effect of 
foreign gas injection on the skin friction and heat transfer of a nearly 
isothermal boundary layer by means of mixing length theory. Numerical 
results are presented for the injection of hydrogen and helium into the 
boundary layer for a temperature of 500° R. It has been found that the 
injection of a given mass of light gas is much more effective than the 
same mass of air in reducing skin friction and heat transfer on the flat 
plate. The same reductions are generally achieved with about 20 percent 
as much hydrogen and about 40 percent as much helium as in the case for 
air. 


INTRODUCTION 


The cooling of aircraft experiencing aerodynamic heating is becoming 
increasingly necessary as the speeds of contemplated aircraft become 
higher and higher. Even with the best of presently available high- 
temperature materials, the steady-state heating of external surfaces of 
aircraft flying at Mach numbers around 10 and higher will often require 
extensive cooling. High-level heating regions such as fuselage tips and 
wing leading edges require cooling at even lower Mach numbers. Cooling 
may even prove effective in transient heating systems, such as in ballis- 
tic missiles, because the heat absorbing ability of fluids used in cool- 
ing systems, such as water, hydrogen, or helium, greatly exceeds the heat 
absorbing ability of solid materials, without phase change, on a pound 
per pound basis. 

Of the various cooling systems available, mass transfer systems in 
which the coolant is ultimately introduced into the boundary layer in 
contact with the aircraft surface appear to be more effective than the 
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more conventional internal cooling systems (ref. l) . The present paper 
is concerned with a transpiration cooling system in which the coolant 
passes through the surface it is protecting before entering the surround- 
ing boundary layer. The boundary layer is considered to be turbulent. 
Analyses and experiments have been performed to determine the effect of 
distributed air transpiration through flat surfaces over which air flows 
in a turbulent boundary layer (refs. 2, 3, and 4 ). These investigations 
revealed that the transpiration process reduces both the skin friction 
and heat transfer associated with the boundary layer. These effects are 
qualitatively the same as the behavior of the laminar boundary layer 
under similar conditions at the surface (refs. 5 and 6). Now it is known 
(refs. 7 and 8) that injection of a light gas into the boundary layer at 
the surface is much more effective than the injection of air in reducing 
the skin friction of a laminar boundary layer on a flat plate. 

It is the purpose of this study to determine if light gas injection 
affects the turbulent boundary layer in an analogous manner. The study 
is confined to the boundary layer on a flat plate and is divided into 
three parts : 

(a) A derivation of the basic fully turbulent boundary- layer equa- 
tions for a binary gas system, considering the diffusion only due to 
concentration gradients, 

(b) An evaluation of modified Reynolds analogies between momentum, 
mass, and heat transfer for a binary gas system, and 

(c) A determination of the effect of foreign gas injection on the 
skin friction and heat-transfer processes of a nearly isothermal turbulent 
boundary layer with the numerical results confined to the injection of 
helium and hydrogen. In this latter part, mixing length theory is 
employed even though it is known to be physically incorrect (ref. 9). 

The pragmatic viewpoint is taken. Mixing length theory has provided use- 
ful approximations in the past (refs. 2, 3, and 10), so it may be expected 
to do so in the present case as well. 


SYMBOLS 


Cj. local skin-friction coefficient 

c-n specific heat at constant pressure per unit mass 

ir 

C constant of integration 

D-^ molecular diffusion coefficient for a binary mixture 

F injection mass flow per unit area divided by the mass flow per unit 

Pw v w 

area of stream just outside the boundary layer, 
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g symbol for exponential term in equation (68) 
h heat-transfer coefficient 

i enthalpy per unit mass 

K mixing length constant 

m symbol defined by equation (65) 

M molecular weight 

n symbol defined by equation (65) 

p pressure 

M* c p 

Pr Prandtl number of mixture, 

Pr* Prandtl number as defined by equation (^3) 


<lc 

R x 

% 

R e 

Sc 

Sc* 


heat-transfer rate per unit area by molecular conduction 


Poo U oo X 

Reynolds number, ~ n 

*00 

^ PM 

Reynolds number, — 77 

MX) 

v p ™ U °° 9 

Reynolds number, — £ 


Schmidt number, 


12 


turbulent Schmidt number, 


e v 

P e d 


St Stanton number, n u c-^ 

Poo^co p^ 

T temperature 

u mass velocity parallel to plate surface 

~ u 

u dimensionless velocity, ^ 

v mass velocity normal to plate surface 

x distance along plate from leading edge 

y distance normal to plate surface 


4 


MCA TN 4149 


a, (3 defined by equations (8l) and (82) 

T defined by equation (44) 

5 boundary-layer thickness per length of run 

eddy diffusion coefficient 

6y eddy viscosity 

6-, eddy thermal conductivity 

A 

F 

? Syi 

0 momentum thickness 

7\ molecular thermal conductivity 

p molecular viscosity 

p density 

a temperature recovery factor 

t frictional shear stress 

lo mass fraction of foreign gas 

Superscripts 

( ) temporal mean value, see equation (8) 

* randomly fluctuating value 

Subscripts 

a condition at interface 

a- condition at interface, on laminar sublayer side 

a+ condition at interface, on turbulent side 

w condition at surface of plate 

oo condition at outer edge of boundary layer 
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l 


light gas 


2 


air 


o 


zero injection condition 


ANALYSIS 


The analysis of this report is quite similar, in principle, to the 
analysis of reference 3* The main difference is that the basic equations 
used in the present analysis apply to a two-component mixture of gases , 
a foreign gas and air, rather than air alone. Air is treated as a single 
gas, having the mean properties of the mixture containing the following 
volume fractions: 0.78 N 2 , 0.21 0 2 , 0.01 A. 

The turbulent boundary layer is considered to occur on a flat plate 
(i.e., dp/dx = 0) oriented parallel to the free stream. It is postulated 
to be composed of two regions: a laminar sublayer, where the transport 

of momentum, energy, and mass is controlled by molecular motion; and an 
outer turbulent region where the transport of properties is controlled 
by eddying motion. At the interface of these two regions, it is required 
that there be a continuity of velocity, temperature, shear, and mass and 
energy flux. The boundary layer is considered to be at steady state. 

Only the diffusion resulting from concentration gradients is included; 
the other diffusion processes are considered to be sufficiently small for 
the purposes of this report so that they can be neglected. 


When the usual boundary-layer order of magnitude argument is 
employed, the laminar transport equations for a bicomponent mixture on a 
flat plate can be derived from the equations for a nonreacting gas mix- 
ture of reference 11 (pp. 698, 498, and 516) . The equations can be 
written as follows: 

Continuity of mass 


Basic Boundary -Layer Equations 



Conservation of momentum in x direction 



( 1 ) 
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Conservation of momentum in y direction 




Conservation of diffusing component 


CjU) 

bx 


^03 S ( -p. Su) 

pu ^ + pv ^ ^ ^pD x2 ^ 


(3) 


(4) 


Conservation of thermal energy 

Si Si S [, 5 t n dto /. _ 1 \ 

pu S + pv 3? " 57 L A 37 p 12 3y ( 1 2> 


+ \l 



(5) 


An equation representing a conservation of total energy, i + —, 
can he obtained by combining equations (2) and (5). The resulting 
equation is 

U 2 \ S ( . S [", St p. So) /. _ . v .... Sul 

^j +pv ^r"j>R + p 12 ( 1 2 ^yJ 

(6) 


pu S 1 1 + 


The corresponding equations for the turbulent portion of the 
boundary layer can be formed if it is assumed that the molecular trans- 
port terms, the right members of the above equations, are equal to zero. 
This implies that the left members represent an instantaneous convection 
of the corresponding properties and that molecular mechanisms are negli- 
gible. The instantaneous values of the dependent variables are then 
expressed as the sum of a term that is invariant with time and a term 
that varies rapidly with time in a random fashion. Thus, for example. 


u = u + u 1 
v = v + v f 
p = p + p* > 

tO = (0 + (0 * 

i = i + i* 


(T) 


The expressions (j) are then substituted into the left members of 
equations (l) , (2) , (4) , and (6) which in turn have been set equal to 
zero. When averages of fluctuating quantities are defined as 
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A' 



(8) 


where r is an interval of time sufficiently large to obtain a meaning- 
ful average, and when the following orders of magnitudes are applied, 


v, p'u 1 , p'v', etc. ~ 0(5) 
p'u'v', p'u'u', etc. ~ 0(S 2 ) 
I, u, p, u, ^ ~ 0(1) 

A „ P*( T ) ~ P*(~ T ) 

^ °( 5 >> 2 t 




(9) 


0, etc, 


( 1 ) ' , p', v ’ , u', i ’ «1 


where 5 « 1 represents the boundary-layer thickness per unit of length 
of run, there result the following turbulent boundary-layer equations 


8x 


(pu) 




+ p’v’) = 0 


(10) 


du / — ; — r\ 8il 

pu _ + (pv + p'v') ^ 


pu 22 + (pv + p'v') |£ 
dx dy 



(11) 


( 12 ) 



(pv + p'v') 



_d_ 

8y 


dT — dio /— — \ , — 

e A^ + pe a ^ v 


du 

8y. 


: ] 


( 13 ) 


where 


_ u»v» 

£v " _P 

v 7 ? 7 I 

€ A - " pc P df/dy f 


v 1 ^* 

Ow/Sy 


( 14 ) 
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Boundary Conditions 

The boundary conditions imposed on the above equations are at y = 0: 


as 


oo : 


u = 

0 

T = 


0) = 


H* 

II 

ii(T w ) 

^ 2 W 

I 2 ( lw ) 

u = 

^00 

T = 

T 

^00 

0) = 

0 

II 

8 

H 

•H 

1 1 ( loo) 

1 2 00 

i-2 ( "^oo ) 

3/fc 

y = 0 (eq. 


( 15 ) 


(16) 


nation. un_i_y m one ui ouiiaoc ^ v ~w 

prescribed, since then the partial pressure of a vapor over a liquid is a 
function of the liquid temperature. If the light gas is injected into 
the boundary layer, the concentration at the surface, Ww, will depend on 
the diffusion rate through the boundary layer as well as the injection 
rate. This dependence can be seen from the following arguments. The 
average mass flow per unit area at the surface can be written as 


Pw v w - (Pi v i)w + ( P2^2 )w 


( 17 ) 


At the steady- state conditions, after air has diffused into the surface 
and established an equilibrium concentration, there is no more net air 
flow through the surface. Therefore (v 2 ) w = 0 and 

Pw v w = (Pi v i)w 

Equation (l8) can be rewritten as 

(Pl W + P 2W )V V = Pl V V % 


(18) 


s 
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since p = p x + p 2 , and dividing through by p^ results in 


K 


jX + (jX}” = (SX 




Since 



by definition, and 


U) 


pi 

p 


(1 “ “w) V W = “w( V l w - V w) 


(19) 


Now v n - v TT defines the diffusion velocity of the injected gas and is 
- L w w 

represented by 



Therefore 


(! - w w) v wPw = -Pw D 12 W (^) W (20) 

Equation (20), then, is a more useful definition for the surface 
concentration, uv, than (15) for a nonevaporative system. 

The above system of equations, together with the boundary conditions, 
is intractable because of its complexity and a lack of knowledge about 
the functional dependence of the eddy coefficients. At present, then, 
approximate methods are necessary. 


Simplified Boundary-Layer Equations 


Experience with the theory for turbulent boundary layers on flat 
plates (e.g., refs. 2 and 3) has shown that for practical purposes, varia- 
tions of dependent variables at a fixed value of x can be found by 
treating the boundary-layer equations as total differential equations in 
y alone, the terms containing differentiation with respect to x being 
neglected. In this process, ultimately, variations in the x direction 
are introduced by the x variation of the shear and heat transfer at the 
wall. This is the technique employed in the present analysis also. The 
simplified boundary-layer equations can be written for the laminar 
sublayer as 
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_d_ 

<3y 


and. 


du _ _d_ / du\ 
d y " dy ^ d y) 


dw _ _d_ ( dtA 
Pw v w dy " dy \_ pD i2 dy/ 


. v d (* , u 2 ^ _ d 

^ w ^ V “/ d7 


P v = Pw v w 


(21) 

du\ 

<3y/ 


(22) 

dw\ 
: ayy 


(23) 

H 

•H 

as 

U + “t* fa] 

(24) 


The corresponding equations for the turbulent portion of the boundary 
layer can be written as 


and 


^ (pv + p*v‘) = 0 or pv + p’v* = const. = p w v w 


„ du _ d [ _ du 
^ w dy dy \ v dy. 


da) _ d 




CLO) 

^ w dy ~ dy 

Pw v w ^ (* + If) = dy [ e A dy + pe d- dy hi " + U€v 

The superscript bars are no longer needed for clarity. 


(25) 


(26) 


(27) 


( 28 ) 


Derivation of Reynolds Analogies 


Before one can solve for the skin friction and heat transfer 
associated with the boundary layer, it is necessary to relate the local 
temperature and concentration with the local velocity. This facilitates 
later integrations and also permits analogies between momentum, mass, 
and heat transfer. 
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Analogy between skin friction and diffusion .- In the sublayer, the 
relation between momentum and mass transfer can be determined from equa- 
tions (22) and (23). If each is integrated with respect to y and the 
boundary conditions at y = 0 are imposed, there results 


Pwv w u + t w = p ^ (29) 

and 

Pw v w( w " !) = P D i2 (30) 

If the right and left members of equation (30) are divided by the 
corresponding members of equation (29) , and the resultant equation 
integrated employing 


Sc 


P D 12 


const. 


for simplicity. 


there results 
m - 1 


r F5 + (c f /2) -j 
|_Fu a + (c f /2)J 


(3D 


(32) 


In equation (32) the arbitrary constant of integration was evaluated at 
the interface and the variables were made dimensionless. 

In the outer turbulent portion, equations (26) and (27) are 
integrated with respect to y and subjected to the requirements that 
the velocity, concentration, shear, and diffusion rate be continuous 
across the interface, namely 

Ua " ■ Ua+ ' O 1 3y) a - ' ( Sv ay) a+ 

“ a - " “ a+ ’ C 0 ” 12 i) a - = ( p6 i t) a+ 

There results 

Pw v w u + T w = e v dy (3^) 


(33) 


Pw v w( a > - 1) = P^d f 


(35) 
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Again, as in the sublayer, dividing equation (35) by equation (34), 
performing the integration with 


Sc T = — — = const. 

p £ d 


(36) 


and applying the boundary conditions at y — > oo yields 


1 - 0 ) 


'Fu + (c f /2)~ 

_F + (c f /2) 


(37) 


When the value of u> a at the interface is eliminated between equa- 
tions (32) and (37) and. the wall condition is evaluated, there results 


1 - ov - 


(c f /2) Sc 

[Fu a + (c f /2)] 

| Sc * - Sc 

1 

[F + (c f /2)] 

| Sc* 


(38) 


which represents an analogy between diffusion and momentum transfer. 

Analogy between skin friction and heat transfer .- The derivation of 
the relationship between heat and momentum transfer is more complex than 
the work of the preceding section because the effect of the simultaneous 
mass transfer must be isolated. When the Prandtl and Schmidt numbers 
(laminar or turbulent) are equal, it is possible to consider total 
enthalpy only; however, when they differ, it is necessary to consider 
both temperature and concentration. The assumption of equal Prandtl and 
Schmidt numbers is made in the turbulent portion of the boundary layer. 

In the laminar sublayer integration of equation (24) together with 
incorporation of equation (30) results in 

Pw v w(^ + t) = A ^ + Pw v w( w - !)( i i - i 2 ) + + C (39) 

Now, the enthalpy term on the left, although largely dependent on T, is 
also dependent on the concentration go. The members on the right contain 
terms which are dependent on temperature alone, dT/dy and i x - i 2 . To 
make the terms of the equation more consistent, the enthalpy term i of 
the left member is expanded as 

i = (oi 1 + (l - co)i 2 


and equation (39) becomes 
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Pw “ i l w + 2) ~ <lc w “ ^ 


dT , du 
dy + ^ dy 


where the arbitrary constant evaluated at y = 0, is 


c = PwVi v + <lc w 


(40) 


(4l) 


The independent variable in equation (4o) is changed from y to u by 
dividing corresponding members of equation (4o) by those of equation ( 29 ). 
In addition, the first term on the right of equation (40) is multiplied 
and divided by Cp_^. There results 


Pw v w [l 1 " ii w + (u 2 /2)] - q Cw 
Pw v w u + T w 


A di i 
pcp^ du + u 


(42) 


Note that the condition of the equality of temperature of the individual 
components and the mixture is employed here. It is also to be noted that 
a peculiar type of Prandtl number appears in equation (42), where the 
normal term Cp is replaced by Cp . If, for mathematical simplicity. 



( 43 ) 


is considered constant across the sublayer, equation (42) can be 
integrated directly by use of an integrating factor. When the boundary 
conditions at the surface are employed, there results 


In " 


u c 


In + 

2 


o Uco 


u[fu + (c f /2)J [fu + (c f /2)J 


(2 - Pr*)F 2 


2-Pr 


(c f /2) 

(2 - Pr*)F s 


(» * 








^ * 
-Pr r 


Fu + 


©r 


= r 




To express equation (W*-) in terms of temperature, it is convenient to 

let c-r, = const. At the sublayer side of the interface, then, 

*1 


T r a - ( S a 2 / 2 ) u oo 2 . T 
T a- - 5T + 


w 


( 45 ) 
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It is now necessary to determine the temperature on the turbulent 
side of the interface, T a+ , which will he matched to T a - of equa- 
tion (^5). By integrating equation (28), there results 

Pv^i + f)=c x § + pc a | (1, - i 2 ) + us v & + C (46) 


Continuity at the interface of heat transfer due to diffusion and 
conduction, of shear, and of enthalpy and velocity requires the constant 
of integration of this equation to be the same as that for equation (4o) , 
namely, equation (4l). To simplify the equation, di/dy is evaluated as 

di _ 8i 
dy " 8T 

By means of the definition of i, i = ijW + i 2 (l - w) , it can be shown 
readily that when the turbulent Prandtl and Schmidt numbers are equal 

= Jl (48) 

e A P £ d 

that substitution of equations (4l) , (47), and (48) into (46) results in 


dT di 
dy dm 


dm 
m <*y 


(47) 


Pw v ) 



fA di 

c p dy 


ue 


du 
v dy 


(49) 


The independent variable is changed from y to u 
equation (34); thus 

Pw v w[i - ^ + ( u2 / 2 )] " % 
p w v w u + r w 


by dividing through by 


di 

du 


+ u 


(50) 


when, for simplicity, the turbulent Prandtl or Schmidt number is set 
equal to unity. When equation (50) is integrated with c p ^ set constant 

and the boundary conditions (l6) are applied, there results 


F + (c f /2) 

F l 

^ c P 2 ^oo “ C P X T W + ( U 00 2 / 2 )J 

- (9c w /P a 

> U J 

Fu a + (c f /2) F 

[ w a c p 1 T a+ + C 1 “ “a) c p 2 T a+ “ c p 1 T w 

+ (u a 2 /2)] 

- (lc w /poo U oo) 


(51) 


When equation (5l) is solved for T a+ and equated to the right member of 
equation (45), there results, after some manipulation. 
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<Jc 


c p (T w - Tj - o(n*/2) 


W 


Poo u o/ 


Fu a +(c f /2) 1 c p x |Fu a + (c f /2) 


F+(c f /2) c p a F+(c f /2) [Fu a + (c f /2)] Pr _1 [f+ (c f /2)J 


(c f /2) Pr 


(52) 


where 


F = 


Pw v w 

Poo^oo 


and 


o = 


F+ (c f /2) " 

Cp a , 

2 Ua 

[F+ (c f /2)‘ 

] 4 

[p%+(c f /2)][p+(c f /2)] 

1_Ua Fu a + (c f /2) 

' C P X 


F 


(2 - Pr*)F 2 


2(c f /2) 2 ~ Pr *[Fu a + (c f /2)] Pr " X [f+ (c f /2)] _ F+ (e f /2) 


(2 - Pr )F 2 


“ u a 


Fu a + (c f /2) 


(53) 


If we let 


!c w = h( T w “ T oo “ a 2^ 


(54) 


and noting that c-q = c-p , then 

*2 *oo 


St 


F/(c f /2) 


r / 2 


F(1 - ua) C P P 


Fu a + (c f /2) c p x 


F+ (c„/2) 


Fu a + (Cf/2) 


[Fu a + (c f /2)J 


(55) 


Pr 


(c f /2) 


Pr 


Equation (55) represents the modified Reynolds analogy between heat trans- 
fer and skin friction when turbulent Prandtl and Schmidt numbers are 
unity* The term cr, defined in equation (53) > is seen from equation (5*0 
to be a temperature recovery factor. 


Effect of Light-Gas Injection on Skin Friction of 
a Near Isothermal, Turbulent Boundary Layer 


In order to evaluate the skin friction exerted by a turbulent 
boundary layer composed of a mixture of a light gas and air, the light 
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gas being transpired at the surface, it is necessary first to integrate 
the momentum equations ( 29 ) and (34) across the boundary layer. Because 
these equations contain the viscosity and density of the mixture 
(through p and e v ) , integration of these equations in general must con- 
sider both the concentration and temperature variation across the boundary 
layer. In the preceding section it was shown that the local concentration 
and temperature could be related to the local velocity; therefore, it is 
possible to relate the local viscosity and density to the local velocity. 
Thus, integration of (29) results in the following equation written in a 
dimensionless form 


, Poo u tt7 = r u (p/O dS 
y ^00 J 0 Fu + (c f /2) 


(56) 


Before equation (34) can be integrated, it is necessary to express 
the eddy viscosity, e v , in terms of a mixing length. The usual expression 
obtained by employing the Prandtl mixing length concept is 


ev 


PKV 


du 

ay 


Employing equation (57) with equation (34) yields 


Pw v w u + T w 



(57) 


(58) 


When the root of both sides of equation ( 58 ) is taken, K is assumed a 
constant, and the resulting equation integrated from u a to u, there is 
obtained in dimensionless form 


R y = R y a ex P 



(p/pj 1/2 du 
«V(c f /2) + Fu 


(59) 


For the present, let it be considered that the boundary layer is nearly 
isothermal so that in integrating equations ( 56 ) and (59) it is not 
necessary to account for temperature variations. The viscosity and 
density thus become functions only of the concentration. 


In the previous sections, the fluid properties within the sublayer, 
such as Schmidt number, have been treated as constants to simplify the 
integrations. The thought behind this simplification is that average 
values of the fluid properties in the sublayer should suffice for the 
calculations of the concentration, momentum, and energy in the sublayer. 
When this idea is carried over to the integration of equation ( 56 ), p/p^ 
is replaced by (p/p ) , and there results 

' 00 cnr 
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or 



In 


Fu + (c f /2) 

c f /2 



/V\ 1 Fu a + (c f /2) 

W av F 1 c f/ 2 


(6o) 


(6i) 


when R y is evaluated at u a , the interface. 

In order to integrate equation (59) in the turbulent portion of the 
boundary layer, it is convenient to employ the density variation 


= T i (62) 

p oo [(Mg/M-J - lju + 1 

This equation applies to an isothermal boundary layer on a flat plate 
and has been obtained from Dalton's law for partial pressures and the 
perfect gas law. Equation (62) must be expressed in terms of the local 
velocity to facilitate integrating equation (59)* From equation (37 ) > 
when Sc * = 1, 


Fu + (c f /2) 

1 - u = , . 7 

F + (c f /2) 

Combining equations (62) and ( 63 ) results in 


(63) 


where 


JL = 1 

Poo m - nu 


(64) 


_ Mg _ /Mg _ \ c f / 2 

" M i " \ M i " / F + (c f /2) 

> 

n _ Ak ^ F 

V «1 7f+ (c f / 2 ) J 


( 65 ) 


Use of equation (64) permits integration of the argument of the exponent 
in equation (59) to yield for n > 0 (light-gas injection) 


R\ 


Zn 


R 




-K 

sl¥n 


-2Fnu - n(c|./2) + mF -2Fnu a - n(c|./2) + mF 

; Sin : — : 

n(c f /2) + mF n(c f /2) + mF 


sin 


(66a) 
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and for n < 0 (heavy-gas injection) 


in y = 

tanh 1 j 

1 -n[Fu + (= f /2)] _ wl 

1 “ n [ FG a + (Cf/2)] 

R y a 4-Fn 

L J 

F(m - nu) J 

F(m - nu a ) 


(66b) 


Equations (6o) and (66) represent the expressions for velocity 
distribution in the boundary layer. 

These velocity distributions will be used to obtain skin-friction 
expressions. Note that only the turbulent velocity distribution is used 
in the integration for the momentum thickness. From reference 3, the 
skin-friction integral equation under conditions of surface injection is 


c f _ d rYpV/! dy de 

t * F = to J qsy u(1 - u) aa du - s 


( 67 ) 


where the integral term represents the momentum thickness, 9. 
From equation (59) it can be deduced that 


dy K (P / Pj 

du 


1/2 


exp 


J( c f /2) + Fu [ \ J( c f /2) + 


K 


X 


u (p/pj l/2 du 


Fu 


( 68 ) 


If the exponential term is represented by the symbol, g, the right member 
of (68) is y a (dg/du). Substitution of equation (68) into the integral 
representing 9 results in 


e ' " a) i d5 


(69) 


When use is made of equation (64) and the independent variable is 
changed to g, equation ( 69 ) becomes 


r*(% 

4(o) 1 


e - y. 




m - nu 


( 70 ) 


When equation (70) is integrated by parts, there results 

,;L (l - 2u) (m - nu) + nu(l - u) 


0 = y« 


u(l - u) 

— g 


m - nu 


o J 0 


g 


(m - nu) 2 


du 


(71) 
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Since g(o) and g(l) are finite, the first term vanishes because of 
the boundary conditions u = 0 at the wall and u = 1 at the outer edge 
of the boundary layer. From (68) and the definition of g 


g du 


J(c f /2) + Fil 

K(P/PJ 1/2 


~ nu J( c f /2) + Fu 


K 


dg 


(72) 


Using (72), equation (71) becomes 
e ■ 


g( 1 ) s /(c f /2) + Fu (m - 2mu + nu 2 )dg 

3/2 


; g(o) 


K(m - nu) 


(73) 


Again, integrating equation (73) "by parts (noting that g(l) » g(0) and 

for c f /2) + F » (c f /2) + F) and retaining only the largest terms, 
yields 


R e 



g(i) It L , 

K n/ 2 


F 


( 7*0 


Equation (74) represents the relationship between the local skin-friction 
coefficient and the local Reynolds number based on momentum thickness. 

To obtain the relationship between the local skin-friction coeffi- 
cient and the Reynolds number based on length of run, it is necessary to 
employ equation (67) in the form 


ff de _ 

2 + 11 ~ dx dR x 


(75) 


On integrating equation (75) there results 

' R 0 dR, 


R> 


-l 


e 


o (c f / 2 ) + F 


if the start of the boundary layer is at R x = 0. By the use of 
equation (74) equation (76) can be transformed to 


R 


x 


r=f/ 2 a{ny a [g(i)/K] J( c f /2 ) + f) 

v/ c =00 (c f /2) + F 


With the use of the relations for R v , g(l), and — r =r, 

y a' d[(c f /2) + Fj 

repeated integration of equation (77) by parts results in 


(76) 


(77) 
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R x = 




Kj(= f /2) 


+ F 


c f / 2 


'f ' ^ smaller 

+ additional 
terms 


( 78 ) 


When the boundary values are imposed, there results 


Ry = 


Ry g(l) 


K 


J( c f/ 2 ) 


( 79 ) 


+ F 


Equation (79) represents the relationship between the local skin-friction 
coefficient and the local Reynolds number based on length of run. In 
both (74) and (79)> g(l) is defined for light-gas injection (n > 0) 


g(l) = exp) - 


K 

<v/Fn 


_ i -2Fn -n(c f /2) +mF 
n(C|,/2) +mF 


sin 


sin 


-2Fnu a - n(c f /2) +mF 
n(c^/2) +mF 


(80) 


and for heavy-gas injection (n < 0) 


g(l) = exp 


2K 

J-¥n 


tanh 


,F±Z 


(c f /2)] -i/-n[FUa + (c f /2)] 

- tainh ' — 


F(m - nu a ) 


RESULTS AND DISCUSSION 


Before the end results of the analysis could be evaluated numerically, 
it was necessary to determine for the gas mixtures considered the molecular 
transport properties, viscosity, thermal conductivity, and diffusion coef- 
ficient, which can be combined to yield Schmidt and Prandtl numbers. The 
viscosity was determined by the method of Wilke (ref. 12); the thermal 
conductivity, by the method of Lindsay and Bromley (ref. 13 ); and the 
diffusion coefficient, by the method of Hirschfelder, Curtiss, and Bird 
(ref. 11). These methods have been recently summarized and compared with 
other available methods by Carlson and Schneider (ref. l4). It can be 
concluded from the latter paper that the methods used in this report are 
in keeping with the present knowledge of the molecular transport mechanism 
in binary gaseous mixtures. 

The results of these calculations are shown in figures 1 through 5 
for both helium-air and hydrogen-air mixtures at a temperature of 500° R. 
Figure 1 shows the variation of the viscosity of the mixture with the mass 
fraction of the light gas. It is observed that the viscosities of the 
mixtures vary continuously as the mixture concentration is changed; however. 
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it is interesting to note that the viscosity of the helium-air mixture 
can exceed the viscosity of the pure individual constituents. The corre- 
sponding thermal conductivity of the gas mixtures is shown in figure 2. 

The results of these figures are combined with the specific heats of the 
mixtures, Cp = Cp^oo + Cp^(l - go) , to form the mixture Prandtl numbers 

shown in figure 3* It should be noted that small concentrations of the 
light gas, 0 ) » 0.1, markedly reduce the Prandtl number of the mixture 
from the value of either of the pure constituents in both mixtures. 

Values as low as 0.42 and 0.46 result for the hydrogen-air and helium- 
air mixtures, respectively. The Prandtl number, Pr*, defined by equa- 
tion (43), in which the specific heat used is that of the light gas, 
rather than that for the mixture, is shown plotted in figure 4. It is 
seen that Pr varies considerably with variations of the mass fraction 
of the light gas, especially at the lower values of go. The importance 
of this large variation of Pr* is discussed later in conjunction with 
the discussion of the heat-transfer results. The Schmidt numbers for the 
gas mixtures are shown in figure 5; the diffusion coefficient used for 
these calculations was determined by means of the equations of reference 11 
for a temperature of 500° R and a pressure of one atmosphere. It is seen 
that the Schmidt number varies considerably with the mass fraction of the 
light gas, the increase being largely due to the decrease in density of 
the mixture. 

It is recalled that in the Analysis, integrations across the laminar 
sublayer were performed considering the Schmidt number, Sc, the Prandtl 
number, Pr*, and the viscosity ratio, constant and equal to the 

average value in the sublayer. This was done to simplify the integrations 
and was thought to be justified by the thinness of the sublayer so that 
only small changes in the transport parameters were expected. To check 
this, curves such as shown in figure 6 were evaluated from the numerical 
results. In this figure there is shown the surface and interface concen- 
tration of the light gases as a function of the injection rate parameter, 

£, which is equal to the local dimensionless injection rate, F, divided 
by the local skin-friction coefficient, c^/2. These curves apply to a 
position where the length Reynolds number is 10 million. It is seen that 
the concentrations of the light gases change about 20 percent through the 
sublayer over most of the £ range shown. From figures 4 and 5, it can 
be seen that a 20-percent change in concentration, depending on the value 
of Gl), can produce large changes in Pr* and Sc. It became necessary to 
investigate what effects these changes could produce in the end results. 

In the determination of go^ from equation ( 32 ), where GJ a is obtained 
from equation (37) > the procedure followed was to estimate the Schmidt 
number, calculate and then determine a new Schmidt number as 

Sc = sc(^^2) 

This Schmidt number was then used in equation ( 32 ) to re-evaluate (%• 

The iteration procedure converged rapidly and required only two steps . 
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A variable Schmidt number which was linear with concentration was used 
to integrate the equations leading to equation (32). The latter results 
were more complicated algebraically and were not significantly different 
from those of the simpler equation (32). The additional complications of 
permitting the Schmidt number to vary, therefore, were not warranted. 

Once the concentration variations across the sublayer were evaluated, 
it was possible to determine the degree of approximation involved in using 
average Pr* and across the sublayer. It was found that using the 

Pr* corresponding to the extreme values of w a and 0% in equation (55) 
produced only about a 3-percent variation in Stanton number, with the 


significant. Also, substituting the extreme value of viscosity associated 
with co a and uv resulted in, at most, a variation of 2 percent in the 
skin-friction coefficient when plotted against length Reynolds number; 
again this is not deemed significant. It can be concluded that the 
simplified mathematical steps of the analysis do not introduce large 
errors and are adequate for an exploratory type of analysis such as this. 

Besides providing information for assessing the accuracy of the 
calculations, the curves in figure 6 provide interesting information 
from the physical viewpoint. The curves of (%■ show how the concentra- 
tion of the light gas at the surface increases with increased injection. 
Only at the highest injection rates shown does approach unity or the 

all light-gas condition. At the low injection rates, £ < 1, the gas at 
the surface is composed mostly of air. 

Numerical evaluation of equation (79) employing equations (6l) 
and (80), with K = 0.392 and u a = 13.1 >s/c f /2 ( ref * 3), results in the 

values of skin-friction coefficient shown in figures 7(a) and 7(b) for 
helium and hydrogen injection into an isothermal boundary layer. The 
parameter of the curves is the term £ and was chosen because the condi- 
tion of a fixed £ results in nearly a constant surface temperature with 
transpiration cooling. The £ = 0 curve is from reference 3* I n general, 
it is noted from either figure that increased injection, larger £, causes 
the skin friction to diminish. This lowering occurs quite uniformly over 
the entire range of Reynolds numbers shown. 

The effect of injection in reducing the skin friction is shown more 
clearly in figures 8(a) and 8(b). Here, the ordinate represents the 
ratio of the skin-friction coefficient to its corresponding value at 
zero injection. The abscissa is the ratio of the local dimensionless 
injection rate divided by the local skin-friction coefficient correspond- 
ing to zero injection. From the numerical results of reference 3 it can 
be shown that these dimensionless groups are convenient coordinates to 
show the effect on skin friction due to air injection in nonisothermal 
boundary layers because the effects of Mach number, wall to free-stream 
temperature, and Reynolds number are largely diminished by use of these 


answers based on 


This would indicate that the effects of the variations in Pr* were not 



roughly halfway between the extremes 
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coordinates. Thus, in an analogous manner, it is considered that the 
isothermal boundary-layer results shown in these figures may be indicative 
of what occurs more generally. Again, the large reductions of skin fric- 
tion by the injection of either of the two gases are evident for each of 
the Reynolds numbers shown. 

A comparison of the relative effectiveness in reducing skin friction 
by injection of helium, hydrogen, and air (ref. 3 ) is shown in figure 9* 
The value R x = 10 7 is picked for this comparison. It is apparent that 
the light gases are much more effective than air in reducing skin fric- 
tion. For instance, a 50-percent reduction in skin friction would 
require mass-flow rates for hydrogen, helium, and air in the proportion 
of 0.21 : 0*39^ • 1«0; or looking at the results in a different way, a 
mass-flow rate of F/(cf/2) Q = 0.80 would cause skin-friction reductions 
of 89 , 62 , and 29 #5 percent for hydrogen, helium, and air, respectively. 


If the restriction that the boundary layer is isothermal is relaxed 
sufficiently to permit some heat transfer without altering to any sig- 
nificant extent the properties of the gas, it is possible to determine 
the effect of injection on the heat transfer by means of the preceding 
calculations of the skin friction. These heat-transfer results can be 
expressed in terns of Stanton number and temperature recovery factor. 

The effect of surface injection of light gases on Stanton number is 
shown in figures 10(a) and 10(b). Here the ratio of Stanton number to 
its vaJ-ue for zero injection is plotted for three Reynolds numbers against 
the dimensionless injection rate, F, divided by the Stanton number corre- 
sponding to zero injection, St Q . This coordinate system, as in the 
coordinate system used for skin friction, has been found convenient for 
air injection to diminish the effects of Mach number, Reynolds number, 
and ratio of wall to free-stream temperature. It is believed that this 
coordinate system may tend to generalize the present results to conditions 
that differ from a nearly isothermal boundary layer. 

It may be seen from figures 10(a) and 10(b) that increasing the 
injection rate from a zero value at first increases the Stanton number 
slightly. At higher injection rates the Stanton number is reduced in a 
manner similar to that of the skin friction. The slight increase in the 
Stanton number for the low injection rate can be explained by examining 
equation (55)* The denominator of the right member is seen to be composed 
of the sum of two terms. The first term, called a and given by 


F(1 ~ u a ) = S(1 ~ u a ) 

Fu a + ( c f/2) £u a + 1 


( 81 ) 


can be considered to be proportional to the resistance to heat flow of 
the outer turbulent portion and, similarly, the second term 0 where 
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' F + (c f /2) “ 


[Fu a+ (c f /2)] Pr 

Cp l 

Fu a + (c f /2) 


. . .Pr* 

( c f/ 2 ) 



( 82 ) 


is proportional to the resistance to heat flow of the laminar sublayer. 

For instance, if u a = 0, and there is no sublayer, (3 has the value zero 
and all the resistance to heat flow occurs in the outer turbulent portion. 
If u a = 1, it is seen that the converse is true, that is, all the resist- 
ance to heat flow would be due to the laminar portion of the boundary 
layer because a is now zero. Now then, what happens at very low injec- 
tion rates? When F, or £, is small, (3 can be rewritten by expanding 
the last bracketed term in a series, thus. 



Pr*5u a 


Pr*(Pr* ~ 1) 


£ 2 G a 2 



(83) 


or 

c p C P 

lim p = Pr *^a = Pr av c^T“ 

£ 0 p i tav 

It is seen that for small injection rates, the resistance of the sublayer 
to heat transfer is proportional to the Prandtl number, evaluated at the 
average concentration in the sublayer and to the ratio c P a / c P av * No ' fce 

that for light gases c p /c p <1 since co a < w av and c p > c p . To 

y a! -^av *1 2. 

illustrate the behavior of a and p numerically, an example is considered 
using equations (8l) and (83)* For hydrogen injection at F/St 0 = 0.046 
and R x = 10 6 , the average concentration of hydrogen, by weight, in the 
sublayer is found from a figure such as figure 6 to be about 3 percent. 

It is noted from figure 3 that this small concentration of hydrogen causes 
a decrease of 30 percent in the Prandtl number of the mixture in the sub- 
layer. For these conditions, the term p is equal to 0.0l6. Now, if 
this drop in Prandtl number did not occur, that is, if the Prandtl number 
and the specific heat remained that of air, then the term p would be 
equal to 0.02^2 and the resistance to heat flow in the sublayer would be 
increased. Since the value of a is 0.0277, the total resistance to 
heat flow through the whole boundary layer is l6 percent less for the 
binary mixture than for the case where the Prandtl number and specific 
heat were the values of air. This tendency to decrease the resistance 
to heat flow, or to increase the Stanton number, appears for the lower 
injection rates to be larger than the normal tendency of the injection 
process alone to reduce Stanton number. At the high injection rates, the 
effect of property changes to increase the Stanton number becomes smaller 
and the effect of the injection process alone becomes larger so as to 
predominate, and causes the decrease of Stanton number with injection 
shown in figure 10 (b). A similar explanation applies for the case of 
helium. 
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A comparison of the relative reduction of heat transfer hy injection 
of hydrogen, helium, and air into the turbulent boundary layer is shown 
in figure 11. The comparison is made for R x = 10 7 . It is evident that 
the light gases when injected into the boundary layer are more effective 
than air in reducing the heat transfer. To reduce the heat transfer to 
50 percent of the zero injection value requires mass-flow rates of hydro- 
gen, helium, and air in the proportion of 0.19 • 0.4 2 : 1.0. As another 
method of comparison, a mass-flow rate of F/St 0 of 0.80 will cause 
heat-transfer reductions of 90, 60 , and 27*5 percent for hydrogen, helium, 
and air, respectively. Included in figure 11 are the experimental data 
of Leadon and Scott, reference 15, for injection of helium and air into 
a turbulent boundary layer for R x = 4xl0 6 , at a Mach number 3.0, and a 
wall to free-stream temperature ratio of 3*30. The curve for air is 
obtained from reference 3 for Mach number zero and wall temperature equal 
to free-stream temperature, that is, a nearly isothermal boundary layer. 
The agreement between theory and experiment is good for air injection. 

If it is assumed that the present theoretical heat-transfer results for 
the helium injection are generalized by the coordinate system of figure 11 
as they are for air, then it is seen that for helium injection the theory 
does not agree closely with experiment but certainly indicates the 
relative advantages of helium injection. 

A comparison of the variation of recovery factor with injection of 
hydrogen, helium, and air is shown in figure 12 for R x = 10 7 . The zero 
injection value, a Q , and the recovery-factor curve for air have been 
obtained from reference 3- The dip in the curves for the low injection 
rates for both hydrogen and helium as contrasted to air may be explained 
by the reduction of the Prandtl number within the sublayer of the boundary 
layer. Comparison with the experimental results of Leadon and Scott, 
reference 15, clearly shows that the theoretical curves do not indicate 
the trends of the experimental recovery factor with injection for either 
helium or air. The disagreement between theory and experiment, evident 
in this figure, emphasizes the need for improving the present theory 
before it can be applied with confidence to predict recovery factors. 


CONCLUDING REMARKS 


A mixing length theory in which it is assumed that both the turbulent 
Schmidt and Prandtl numbers are unity has been developed for the case of 
the injection of foreign gases into the air of a turbulent boundary layer. 
It has been found that: 

1. From the viewpoint of reducing skin friction by injection, helium 
and hydrogen are generally 2.5 and 5«0 times as effective as air on a mass 
basis . 


2. From the viewpoint of reducing Stanton number by injection, 
helium and hydrogen are generally 2.4 and 5*2 times as effective as air 
on a mass basis . 
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3. There are no major effects on temperature-recovery factor due 
to the injection of helium and hydrogen. 

At present, no data exist for evaluating the results of item 1. For 
item 2, however, experimental results generally confirm the high effec- 
tiveness of the light gases in reducing the Stanton number, although there 
does not exist exact quantitative agreement between theory and experiment. 
For item 3, the theory does not even indicate the general trends of the 
experimental temperature-recovery factors. 

Even though mixing length theory is known to be basically incorrect, 
many useful results have been obtained through its use. It has been shown 
in reference l6 in which a mixing length theory was used that taking 
account of variations in the shear distribution across the boundary layer 
and permitting deviations of the turbulent Prandtl number from unity have 
improved the representation of the Stanton number and particularly the 
temperature -re co very factor of a single-constituent gas. From this it 
is indicated that the present analysis should be extended to incorporate 
variations in shear across the boundary layer and deviations from unity 
of the turbulent Schmidt and Prandtl numbers. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Nov. 18, 1957 
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Figure 1.- Viscosity of binary gas mixtures; T = 500° K; method of Wilke, 

reference 12. 
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Figure 2.- Thermal conductivity of binary gas mixtures; T = 500° R; 
method of Lindsay and Bromley, reference 13 . 
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Figure 3 »“ Prandtl number of binary gas mixtures; T = 500° R. 
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Figure 4.- Variation of Pr* with light-gas mass fraction; T = 500° R. 
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Figure 5* “ Schmidt number of binary gas mixtures; T = 500° R. 
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Figure 6.- Wall and interface light- gas concentration variation with 

dimensionless injection parameter. 
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(a) He-Air 


Figure 7«“ Effect of distributed injection on local skin-friction coefficient 
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(b) H 2 -Air 

Figure 7*“ Concluded. 
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(a) He-Air 

Figure 8.- Effect of distributed injection on skin friction. 
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(b) H 2 -Air 

Figure 8.- Concluded. 
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Figure 9*~ Comparison of the effect of distributed injection of He, H 2 , and Air on skin friction. 
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(a) He-Air 

Figure 10.- Effect of distributed injection on heat transfer. 
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(b) H 2 -Air 

Figure 10.- Concluded. 
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Figure 11.- Effect of distributed injection on heat transfer and comparison with experiment. 
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Figure 12.- Effect of distributed injection on recovery factor and comparison with experiment. 
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